Nucleation and Collapse of the Superconducting Phase in Type-I Superconducting 

Films. 
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The phase transition between the intermediate and normal states in type-I superconducting films 
is investigated using magneto-optical imaging. Magnetic hysteresis with different transition fields 
for collapse and nucleation of superconducting domains is found. This is accompanied by topological 
hysteresis characterized by the collapse of circular domains and the appearance of lamellar domains. 
Magnetic hysteresis is shown to arise from supercooled and superheated states. Domain-shape 
instability resulting from long-range magnetic interaction accounts well for topological hysteresis. 
Connection with similar effects in systems with long-range magnetic interactions is emphasized. 



PACS numbers: 74.25.Ha, 05.65.+b, 75.70.Kw 



When submitted to a magnetic field, a type-I supercon- 
ductor undergoes a first-order phase transition between 
the superconducting (SC) and the normal-state (NS) ho- 
mogeneous phases (HP). In the case of films in a perpen- 
dicular field, the transition proceeds through the onset of 
a modulated phase (MP), the so-called intermediate state 
(IS), which consists of an intricate pattern of SC and NS 
domains 0, 0, 0, 0- Such a transition is encountered 
in a variety of quasi-two dimensional systems including 
ferromagnetic thin films with strong uniaxial anisotropy 
[^, |ij , Langmuirpolarized monomolecular layers at air- 
water interface jj, magnetic fluids in Hele-Shaw cells |{|. 
MPs arise from the competition between short-range in- 
teractions associated with positive interface energy and 
long-range magnetic or dielectric interactions. Although 
these interactions have been recognized as a major ingre- 
dient in the description of the dynamics of pattern for- 
mation [J, |tl UiA |ll| , they were not taken into account in 
former studies of interface motion [L^ll3| . An important 
issue concerning MP systems is the role of long-range in- 
teractions in the nucleation process of one phase into the 
other. In a closely related field, the ion-induced nucle- 
ation of the liquid phase in the gas phase of a polar fluid, 
a problem which dates back to the invention of Wilson's 
cloud chamber, this question is still under active debate 
0. 

In type-I SC films, it is well known that the SC-NS 
transition occurs at a magnetic field smaller than the 
bulk thermodynamical critical field H c . This is due to 
the SC-NS interface energy and the magnetic stray field 
energy of the NS domains. The transition field was esti- 
mated using an approximate expression of this magnetic 
energy [l5| . A more accurate prediction of the transi- 
tion field should be obtained in the framework of the 
recently developed current-loop (CL) models 0, 0, 0] 
that made possible the calculation of the magnetic en- 



ergyfor various domain patterns. However, in SC films 
[l|, [l7j], as well as in magnetic systems |6|, a hystere- 
sis loop opens up very close to the MP-HP boundary. 
Two distinct transition fields are found for the appear- 
ance and collapse of domains. Surprinsingly, the origin of 
this magnetic hysteresis still remains an open question. 
Does it arise from the existence of supercooled (SCL) 
and superheated (SH) metastable states? What is pre- 
cisely the role of pinning centers and defects? Metastable 
states were clearly identified in dispersions of micron size 
SC spheres where the small volume reduces the proba- 
bility of heterogeneous nucleation at defects 0, 0] . On 
the opposite, in large size systems like films, SCL and 
SH states are not expected to be observed |2j. In addi- 
tion to magnetic hysteresis, domain-shape hysteresis is 
found: domain shape and pattern are different for rising 
and decreasing field. The interplay between this topolog- 
ical hysteresis and the magnetic hysteresis at the MP-HP 
boundary is not well understood. Valuable insights into 
this question are expected from the study of the stabil- 
ity range of different domain shapes, a task that recent 
models have made possible 0, . 

In this Letter, we discuss the physical origin of mag- 
netic and topological hysteresis close to the MP-HP 
boundary. The two transition fields are shown to cor- 
respond to the rupture of metastable states. They are 
used to determine the critical radius for the nucleation 
and collapse of the SC phase. We show that topological 
hysteresis very likely originates from domain-shape insta- 
bility arising from long-range interactions. The theoreti- 
cal analysis of metastable states and domain-shape insta- 
bilities is carried out in the framework of the constrained 
current-loop (CCL) model which was successfully devel- 
oped to take into account screening by superconducting 
currents in SC films p|. 

The IS domain pattern in SC films was studied with 



FIG. 1: Intermediate state pattern in the 10 /im thick film 
close to the transition to the normal state. The superconduct- 
ing (SC) domains appear in white and the normal phase in 
gray, (a) In rising magnetic field, SC bubbles and lamellae are 
observed {H/H C =Q.TJ at T=1.92 K with H C (T = 0)=282 G), 
(b) In decreasing field, only lamellar domains are observed 
(H/H c =0.68, T=1.85 K). 



the high-resolution Faraday microscopy imaging tech- 
nique 20]. The samples consisted of indium films with 
thicknesses 0.6, 1.1, 1.5, 2.2, 10.0 ± 0.1 and 33 ± 3 
pro.. They were placed in an immersion-type cryostat in 
pumped liquid helium. The optical setup is similar to 
a reflection polarizing microscope 0. The samples were 
zero-field cooled then subjected to a perpendicular mag- 
netic field. 

The SC-NS phase transition is found to be hysteretic. 
Distinct transition fields are found in rising and decreas- 
ing applied field. They are associated with different mor- 
phologies of the SC domains. Figure 2] shows typical 
IS patterns in the 10 /Ltm thick film near the MP-HP 
boundary. In rising field lamellar and bubble-shape SC 
domains are observed (Fig. da)). The length of SC 
lamellae decreases until they are reduced to SC bubbles. 
The transition to the homogeneous NS phase results from 
the collapse of these bubbles whose diameter is 6-7 /jm. 
In contrast, with decreasing field, the lamellar pattern 
appears abruptly in a very narrow range of field (Fig. 
IHb)). The associated magnetic hysteresis is displayed 
in Fig. The area fraction of the NS phase p n , de- 
termined from magneto-optical images, is plotted as a 
function of the reduced applied field h — H/H c . The 
magnetic field H n in the NS domains is related to p n 
by flux conservation p n H n — H. There is a clear devia- 
tion from p n = H/H c shown by the dashed line, which 
means that H n < H c . The transition to the NS (p n =l) 
is thereby completed at a lower field than H c . However, 
the transition is not characterized by a unique transition 
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hup {hdown) are the 



fields at which the SC domains collapse (appear) when h 
is increased (decreased). 

This magnetic hysteresis was observed in all the stud- 
ied samples. In Fig. h up and hdown are reported as a 
function of the magnetic Bond number Bm = d/2irA(T) 
[l6j |. d is the sample thickness and A(T) is the interface 
wall parameter [2lj. In order to determine whether h up 
and hdown are related to superheating and supercooling, 



FIG. 2: Area fraction of the normal phase p n as a function of 
the reduced applied magnetic field h for the 10 pm thick sam- 
ple at T=1.88 ± 0.04 K (Bm=3.2). Black squares: increasing 
field, empty squares: decreasing field, dashed line: p n =H/H c . 
The experimental nucleation (hdown) and collapse (h up ) fields 
are indicated by arrows below and above the calculated crit- 
ical field h cr it, respectively. 
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FIG. 3: Experimental nucleation field hdown (black squares) 
and collapse field h up (empty squares) as a function of the 
magnetic Bond number. The full and dashed curves are the 
theoretical critical and collapse fields, respectively. 



they were compared to the transition field deduced from 
the free energy of the system. In the framework of the 
CCL model, we calculated the free energy associated with 
the formation of an isolated SC cylindrical domain 
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where asN = AH 2 /8tt is the interfacial tension between 
the SC and NS phases, p = 2R/d is the reduced bub- 
ble diameter. The first term in Eq. ^ is the interface 
energy. The second term contains the bulk magnetic en- 
ergy and the condensation energy. It is negative since 
h < 1. The third term represents the interaction en- 
ergy of the screening current circulating within the bub- 
ble wall. The bubble energy is plotted in Fig. 0] as a 
function of p. The set of curves obtained for different ap- 
plied fields presents the typical behavior of a metastable 
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system. The critical field h crit = —x 1 / 2 + (1 + x) 1 ^ 2 , 
x = 8/ (3n 2 Brnj, is the field at which the free energies 
with and without a SC bubble are equal. An energy 
barrier impedes the nucleation or the collapse of a bub- 
ble at h = h crit ■ Starting from the NS phase (p = 0) 
and decreasing H, the system may stay in a metastable 
state. Nucleation of a SC bubble occurs if p > p nuc i with 



barrier. This suggests a common nature for nucleation 
and pinning centers. 



Pnud = y(l-(l-8/z) 1/2 ), y = n(l-h 2 )/(4h 2 ), 

z = 7r 2 (l — ft, 2 ) 2 Bm/h 2 ; h stands here for the nucleation 
field. The expansion of the SC phase is however limited 
by the long-range Biot-Savart interaction of the screen- 
ing current (p 3 term in Eq. ^| leading to an equilibrium 

bubble diameter po = y (l + (1 — 8/z) 1 ^ . In the fol- 
lowing we assume that the first step of nucleation yields 
circular domains. Their evolution towards the laminar 
shape will be discussed later. 

Starting from the IS state with SC bubbles and rising 
the field, the system may remain in a metastable state 
above h cr it, up to the collapse field corresponding to the 
disappearance of the energy barrier: h co ii = —w 1 ! 2 + 
(1 + w) 1 ^ 2 with w = 2 1 '(it 2 Bm) . The corresponding col- 
lapse diameter is p c 



(l + (l + l/ W ) 1/2 )/(^m). 



h cr n and h co u are compared to hd OW n and h up in Fig. 
CD hdown and h up -values lie below and above h cr it , respec- 
tively. This is consistent with the existence of barriers for 
nucleation and collapse. SCL and SH states are indeed 
observed. Since the hdown fields remain much larger than 
the spinodal limit H C 2/H C « 0.12 for indium, the onset 
of the SC phase proceeds through heterogeneous nucle- 
ation. Thermally activated nucleation can be ruled out 
since the barrier is larger than the thermal energy kT by 
many orders of magnitude. Defects likely act as poten- 
tial wells that locally cancel the nucleation barrier when 
its amplitude is sufficiently lowered by the applied field. 
The nucleation radius R n uci — Pnucid/2 is obtained from 
the nucleation field hdown and plotted in Fig. El^a) in 
units of A. Rnuci is of the order of 2 A showing a varia- 
tion by only a factor 1.5 while Bm is changed by a factor 
30. The nucleation volume is thereby ~ 7rA 2 d. This is 
quite reasonable since A is of the order of the coherence 
length £, which is the typical dimension of perturbation 
of the order parameter when a SC domain nucleates. A 
more accurate description of nucleation should take into 
account the spatial variation of the order parameter at 
the SC-NS interface, but this is beyond the scope of the 
CCL model. 

Considering the collapse of SC bubbles, let us note that 
the ft. up -values lie even above h co ii- As the movement of 
SC bubbles is frozen close to h up , this shift of the SC- 
NS transition beyond the collapse field likely originates 
from the existence of pinning centers that form local po- 
tential wells. For all samples studied they are found to 
decrease the energy of the system by almost the same 
quantity as the potential wells that cancel the nucleation 




FIG. 4: Energy of a single SC bubble in the NS matrix in 
units of 2irosNd 2 as a function of the reduced diameter. The 
_Bm-value of the 10 /im thick film, Bm=3.2, is used. 
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FIG. 5: (a) Nucleation radius (in units of A) obtained 
from the nucleation barrier at hdown- (b) calculated critical 
(dashed curve) and collapse (full curve) diameters, collapse 
diameter obtained from h up (black squares), measured SC 
bubble diameter close to h up (empty circles). 



Let us examine now whether the CCL model, which 
describes well the magnetic hysteresis, also provides a 
good agreement for domain sizes. In Fig. EJb) the av- 
erage diameter of SC bubbles measured close to h up is 
compared to the calculated diameters p cr it and p co ii and 
to the diameter p up . p up corresponds to the p- value for 
which d 2 F/d 2 p—0 at h = h up . The measured diameters 
are in quite good agreement with p co ii and p up , thereby 
indicating that the CCL model accurately describes do- 
main sizes. 

Let us now address the question of domain-shape hys- 
teresis: why does nucleation of the SC phase give rise 
to the lamellar pattern even though the ground state of 
the system close to the critical field is the bubble phase? 
It was suggested that the ramification of the SC phase 
propagating in the NS phase originates from dynamical 
instabilities driven by magnetic field diffusion 0, Il3j |. 
However the role of long-range interactions in branch- 
ing instabilities was later emphasized 0, ^| . We con- 
sider here only the onset of domain formation. We show 
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FIG. 6: Equilibrium diameter of SC nucleated bubble (solid 
line) and critical diameter for elliptical instability (dashed 
line), (a) indium (Bm=3.2);(b) tin (Bm— 387), experimental 
data from tf c =305.5 G A=0.23 fim 



growth of the lamellar pattern in order to reach the equi- 
librium state corresponding to the applied field. 

In conclusion, the SC-NS phase transition in type-I 
SC films exhibits magnetic hysteresis and domain-shape 
hysteresis, which are shown to arise from different phys- 
ical phenomena. Magnetic hysteresis, characterized by 
different values of the collapse and nucleation fields of 
SC domains, is found to be the signature of metastable 
states. Domain-shape hysteresis manifests itself as the 
collapse of bubble domains and nucleation of lamellar 
domains. Bubble-shape elliptical instability provides a 
very likely explanation for this topological hysteresis for 
a broad range of values of the magnetic Bond number. 
An analysis along the same lines would be of valuable in- 
terest for other systems with long-range interactions that 
exhibit similar domain patterns and hysteretic behavior. 



that shape instability arising from long-range magnetic 
interactions very likely accounts for topological hystere- 
sis. From linear stability analysis |2^| the critical diam- 
eter for the bubble elliptical instability is obtained as 
Pinst — 3/h\/2Bm. If the bubble diameter after nucle- 
ation po is larger than pi ns t a bubble evolves into an 
elongated-shape domain, po and pinst are plotted as a 
function of the field in Fig. EJa) for the 10 /xm indium 
sample (Bm—3.2). At the nucleation field shown by the 
vertical bar the nucleated bubble diameter is very close 
to the instability limit, being smaller by only 10 %. 

Moreover our theoretical predictions also provide ex- 
cellent agreement with experimental data obtained by 
muon spin rotation spectroscopy on white tin as 
shown in Fig. Elb). The magnetic Bond number Bm = 
387 is much larger than for indium. The field of disap- 
pearance of the bubble phase (h up ) is slightly larger than 
the calculated collapse field h co u. The crossing between 
Po and pi ns t coincides with the field htrans at which a 
transition from the lamellar to the bubble phase is ob- 
served in rising field. In decreasing field, nucleation oc- 
curs at hdown, below the crossing point. Nucleated bub- 
bles are unstable and therefore the lamellar phase is ob- 
served. 

We propose the following description of the hysteretic 
SC-NS transition. In rising field, the SC lamellar phase 
evolves towards the bubble phase which can remain in a 
metastable state above the critical field. The complete 
transition to the normal phase occurs with the collapse of 
finite-size bubbles at the collapse field or slightly above 
if SC domains are trapped in local potential wells. In de- 
creasing field, the NS phase is a metastable state. Nucle- 
ation occurs below the critical field due to the existence 
of a nucleation barrier. Depending on the magnetic Bond 
number and the value of the nucleation field, nucleation 
may yield unstable bubbles with respect to elliptical de- 
formation. They evolve into lamellae with subsequent 
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